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Bell's theorem implies that the outcomes of local measurements on two maximally entangled sys- 
tems cannot be simulated without classical communication between the parties. The communication 
cost is finite for n Bell states, but it grows exponentially in n. Three simple protocols are presented 
that provide approximate simulations for low-dimensional entangled systems and require a linearly 
growing amount of communication. We have tested them by performing some simulations for a fam- 
ily of measurements. The maximal error is less than 1% in three dimensions and grows sublinearly 
with the number of entangled bits in the range numerically tested. One protocol is the multidi- 
mensional generalization of the exact Toner-Bacon [Phys. Rev. Lett. 91, 187904 (2003)] model 
for a single Bell state. The other two protocols are generalizations of an alternative exact model, 
which we derive from the Kochen-Specker [J. Math. Mech. 17, 59 (1967)] scheme for simulating 
single-qubit measurements. These protocols can give some indication for finding optimal one-way 
communication protocols that classically simulate entanglement and quantum channels. Further- 
more they can be useful for deciding if a quantum communication protocol provides an advantage 
on classical protocols. 



I. INTRODUCTION 



For some tasks in information processing, quantum 
communication channels were proved to be much more 
powerful than classical channels. While Holevo's theo- 
rem [l| states that n qubits cannot encode a message 
of more than n classical bits, quantum channels reveal 
their real power when separate devices need to exchange 
information to jointly perform a task whose result de- 
pends on the data held by any single party. In dis- 
tributed computing, the communication complexity of a 
problem is the cost in communication of the most effi- 
cient solution for that problem [2|. Similarly, the quan- 
tum communication complexity is the minimal amount 
of quantum communication required to accomplish a dis- 
tributed computation. Quantum protocols for some com- 
munication problems, such as the Raz's problem Q] or 
the hidden matching (HM) problem [4], made clear that 
in some cases quantum channels can be exponentially 
more powerful than classical channels. Indeed, the HM 
problem exhibits an n (qubits) versus 2^^^^ (bits) gap 
between the quantum and classical communication com- 
plexity for bounded-error protocols. In the case of the 
Deutsch-Jozsa problem and errorless protocols, the gap 
is even stronger: n qubits versus ^(2") bits [5|. A review 
on quantum communication complexity can be found in 
Ref. ^]. 

A protocol of communication complexity is supposed 
to produce a correct result at least with high probability. 
Let us consider the scenario introduced by Yao 12]- Two 
parties, say Alice and Bob, get a fixed part of the input 
data. Alice's input is an element mA of a finite set Ma 
and Bob's input is an element mg of a finite set Mg. The 
purpose is computing a function /(rriA, jtib) with a prob- 
ability of success close to 1 by some amount of communi- 
cation between the parties. In a broader scenario one can 
just require that the outcome is generated according to a 
given probability distribution depending on tua and ms ■ 



The problem of classically simulating a quantum com- 
munication channel fits into this extended framework. A 
quantum state is prepared by Alice through a procedure 
rriA and sent to Bob through a quantum channel. He 
eventually performs a measurement through a procedure 
mg. The measurement outcome, say /3, is generated with 
a probability, p{P\mA,mB), that depends on both mA 
and ms ■ The task of a classical simulation is reproduc- 
ing the distribution p by replacing the quantum commu- 
nication with a classical communication. This task is a 
particular problem of communication complexity in the 
outlined generalized sense. We call the minimal amount 
of classical communication needed in the simulation clas- 
sical communication complexity of the quantum channel 
(the classical communication is not supposed to be neces- 
sarily one-way). Of course, such a simulation requires at 
least the same communication resources that are neces- 
sary for classically simulating a quantum communication 
protocol with an (almost) deterministic outcome. Thus, 
the previously mentioned results in quantum communi- 
cation complexity sharpen the conceptual differences be- 
tween quantum and classical physics. Indeed, they im- 
ply that any classical description of n qubits requires a 
quantity of resources growing at least as 0(2") [^ (or 
20(v^) jf g^ bounded error is admitted), a property that 
hardly fits into the framework of classical physics, where 
in general the quantity of resources scales linearly with 
the amount of information virtually accessible in an ex- 
periment. Thus, a peculiar feature of quantum physics 
is the necessity of a huge quantity of information that 
is almost completely concealed from a direct experimen- 
tal observation, but it is nevertheless fundamental in the 
description of systems. Indeed the whole information en- 
coded in the quantum state cannot be used for carrying 
information [1], but most of that is fundamental in the 
description of some quantum protocols of communication 
complexity. 

There is an open question concerning the classical com- 



munication complexity of a quantum channel. On the 
one hand, the best known protocol for simulating a quan- 
tum channel uses an amount of resources that scales as 
n2^ [y, Q, even with a bounded error. On the other 
hand, the HM problem gives the lower bound 2^^^' 
for the minimal amount of communication in the case of 
bounded error. At present, no other constraint is known; 
thus one could hope to find a better bounded-error sim- 
ulation of a quantum channel with communication com- 
plexity scaling as 2^. 

As established by Bell's theorem [1(|, correlations of 
outcomes produced by local measurements on entangled 
systems cannot be explained classically without post- 
measurement communication between the parties. The 
problem of quantifying the classical communication com- 
plexity of a quantum channel of n qubits is essentially 
equivalent to the problem of finding the minimal amount 
of communication needed to classically simulate the out- 
comes of local measurements on n Bell states. Indeed 
classical protocols for simulating quantum channels can 
be converted into classical models of entanglement with- 
out affecting the cost of communication. Conversely, a 
classical model of entanglement can be converted into a 
model of quantum channels with a little more communi- 
cation, as shown in Ref. [ll| in the case of a single Bell 
state. A more general proof will be given in Sec. IIIIBl 
where we will show that an increase of communication 
by n bits on average is sufficient for the conversion. 

In this paper, we present three approximate classical 
protocols for simulating bipartite entanglement that need 
a one-way communication equal to the number of entan- 
gled bits (ebits). One of these is a multidimensional gen- 
eralization of an exact model for a single ebit, reported 
by Toner and Bacon in Ref. [I4I . The other two protocols 
are generalizations of an alternative exact model, which 
will be derived here. The accuracy of the protocols was 
numerically tested for a family of measurements and for 
a dimension of the Hilbert space (of each party) between 
2 (one ebit) and 32 (5 ebits). The maximum error is 
less than 1% in the three-dimensional case and increases 
sublinearly with the number of ebits in the tested range 
(Note that the error cannot be bounded, since, as previ- 
ously said, a protocol with bounded error needs at least 
20(vn) i^j^g Qf communication). Our models can be use- 
ful for two reasons. First, they can give an indication 
for finding an optimal algorithm with bounded error, as 
discussed in the concluding remarks of the paper. Sec- 
ond, even if the maximal discrepancy increases with the 
number of qubits, these models or some modified version 
can give accurate results for measurements and states 
involved in some protocols of quantum communication 
complexity. In such a case they would give a proof that 
these quantum protocols do not provide any advantage 
on classical protocols. 

In Sec. ini we review the Toner-Bacon protocol for a 
single Bell state and derive the alternative exact protocol 
starting from the Kochen-Specker hidden variable model 
of a qubit 13] . The three generalized protocols for higher 



dimensions of the Hilbert space and their Monte Carlo 
simulations are presented in Sec. IIII Al In Sec. IIIIBl a 
general method for converting an entanglement model 
into a quantum channel model is presented. Finally, we 
draw the conclusions and perspectives in the last section. 



II. SIMULATION OF BELL STATES WITH ONE 
BIT OF COMMUNICATION 

Two qubits are prepared in the Bell state 



|vI/,) = -L(|l)|2)-|2)|l)) 



(1) 



and each one is sent to two parties, Alice and Bob, who 
perform a local measurement on their own qubit. Bell's 
theorem [lO[ implies that the outcome probabilities can- 
not be explained by a local classical theory and any exact 
simulation of this scenario needs some communication 
between the parties. How much information has to be 
exchanged? Trivially Bob could send the whole informa- 
tion about the measurement he performed, which is ac- 
tually infinite. However, Brassard et al. Q showed that 
a finite amount of communication is sufficient for exactly 
reproducing the outcome probabilities. They presented 
a model that simulates Bell correlations and requires ex- 
actly 8 bits of communication. Steiner [lj| reported a 
different model, which requires 2.97 bits on average, but 
the amount of communication for each particular realiza- 
tion is unbounded. This result was improved in Ref. jll| , 
where the average information was lowered to 2.19 bits. 
Later, Toner and Bacon showed that just one bit of com- 
munication can account for Bell correlations [l^l and two 
bits are sufficient for simulating teleportation of a qubit. 
In the next subsection we review this model. In Sec. IIIBI 
an alternative exact model is derived whose direct gen- 
eralization to N dimensions is one of the three protocols 
described in Sec. IIIIAI 



A. Toner-Bacon model 

Let us consider the previously described scenario with 
two qubits in the Bell state ([Ij. Alice and Bob per- 
form two projective measurements. Bob's measurement 
projects the quantum state into one of two mutually or- 
thogonal vectors of the two-dimensional Hilbert space. 
Let us represent them by the Bloch vectors 61 and 
62 = —61. Alice performs a projective measurement on 
the states ai and 02 = —di. The joint probability of 
having outcomes da and bp is 



P(a,/3) = i(l-a„-6^ 



(2) 



It is possible to exactly reproduce this statistics through 
a classical protocol that uses just one bit of classical 
communication. The protocol is as follows. Alice and 



Bob share two random unit vectors Ai and A2. They 
are uncorrelated and uniformly distributed on the unit 
sphere. Bob generates the outcome /3 such that the vec- 
tor bp € {61, 62} is closest to Ai, that is, 



sgn(6^ • Ai) > 0. 
Bob sends one bit nG{— l,l}to Alice, where 

n = sgn(6^ • A2). 
Alice generates the outcome a such that 



sgn da ■ (Ai +nA2) 



<0. 



(3) 
(4) 
(5) 



This protocol produces the events according to the quan- 
tum probability P{a, f3), as proved in Ref. [12[. Note that 
the effect of n is to change the sign of A2, that is, Alice 
receives the instruction "if n is negative, flip the vector 
A2 to the opposite direction" . 

It is useful to express the protocol in a form that is 
trivial to generalize to higher dimensions of the Hilbert 
space. For this purpose we replace the Bloch vectors with 
vectors in the two-dimensional Hilbert space. While the 
state ([T]) simplifies the notation with Bloch vectors, for 
a generalization to higher dimensions it is better to use 
the state 



I*) = ^(|1>|1> 



|2)|2)) 



(6) 



The two states (|1I6P differ by the local unitary transfor- 

- TT " (2) 

mation e^^"y , performed on the second qubit. 

The projective measurements are represented by two 
sets of orthogonal vectors. 



{\^{)M2)} = Ma, { 1^1), 1^-2)} = A/£ 



(7) 



Alice and Bob measure Ma and Ms, respectively. Each 
vector in the set is associated with one outcome. In the 
simulation, the shared vector Ai is replaced by a vector 
|a;) in the two-dimensional Hilbert space, 



Ai -^ \x) 



(8) 



so that Ai is the Bloch vector of \x). The vector A2 is 
replaced by a set of orthogonal vectors in the Hilbert 
space. 



A2^{|yi),l2/2)}, 



(9) 



so that A2 and — A2 are the Bloch vectors of |yi) and |y2), 
respectively. In this new frame, the Toner-Bacon proto- 
col is as follows. According to Eq. (jSj, Bob generates the 
outcome IV'^) € Mb that is closest to \x) (he maximizes 



l(^/3|a 



Then he sends Alice the index n G {1,2} such 



that \yn) G {|yi), I2/2)} is the vector closest to 1?/;^) [from 
Eq. ([U] . Let us introduce the vector 



|0*)^(0fe|l)|l) + (0fe|2)|2) 



(10) 



for fc = 1,2. It is obtained from |(/)fe) by the complex 
conjugation of the coefficients in the basis {|1), |2)}. Alice 
generates the outcome \4>a) G Ma that maximizes the 
function 

{^l\{\^){A + \yn){yn\}Wa). 

Note that Eq. (O would ask to minimize 
{4ia\{\x){x\ + \yn){Vn\}\(i>a) in the case of state ^. 
It is easy to show that for state ([6]) "minimization" 
is replaced by "maximization" and \4>a) by | </<„). The 
outcomes are generated according to Born's rule, that is. 



P(a,/3|M^,Af£ 



1 



\{<t>o.m? 



(11) 



This new reformulation of the Toner-Bacon protocol 
suggests a very simple generalization to higher dimen- 
sions of the Hilbert space, as discussed in Sec. IIII AI 



B. Alternative model 

An alternative exact classical model for a Bell state can 
be obtained from the hidden variable model of a qubit 
introduced by Kochen and Specker (KS)pJi]. This can 
be seen as a classical model of a quantum channel where 
the communicated classical information is infinite and 
encoded in a three-dimensional unit vector. It can be 
transformed into a model with finite communication cost, 
as shown below. 



1. Quantum channel 

Let us introduce the KS model. It provides a classi- 
cal simulation of the following scenario. Bob prepares a 
qubit in a quantum state represented by a Bloch vector 
h. He sends the qubit to Alice, who performs a projective 
measurement Ma — {0.1, 0,2} with outcome states Si and 
3,2 = — fli. The classical simulation is as follows. Given 
the Bloch vector b, Bob generates a unit vector xi with 
probability 



p{xi\b) = -b-xi9{b- xi), 

TT 



(12) 



where 9{z) is the Heaviside function, which is equal to 
1 if z > and zero otherwise. He sends xi to Alice. 
Given the vector pair Ma = {oi, 02}, Alice generates the 
outcome «„ according to the conditional probability 



P{a\xi;MA) = 0{da ■ fi); 



(13) 



in other words, she generates deterministically the vec- 
tor closest to Xi. This model gives the correct quantum 
probability for the outcomes, that is. 



/ d^xiP{a\xi;MA)p{xi\b) = - (l + a,, • &) . 



(14) 



Let us define the vectors 



2. Entanglement 



2/1 

m 



Xl 
Xl 



X2, 
X2. 



(15) 



It is easy to show that the probability distribution P^ 
can be obtained from the uniform distribution 

p{x,,x2\b) = ^e{b-yi)9{b-y2) (16) 

by integrating out the unit vector X2. The integration 
is trivial in spherical coordinates and is left as an exer- 
cise. This uniform probability distribution and the con- 
ditional probability (J13p still give the correct quantum 
predictions. More generally we get an exact simulation 
of a qubit with the probability distributions 

p{xi,X2\ni,n2;b) = —^9 (nib-yi) 9 \n2b- y2J (17) 

P(a\xi,X2,ni,n2\MA) = 6* [&„ • {niyi +^21/2)] (18) 

for any rii — ±1. For example, with ni — \ and 71,2 ~ 
-1, Eqs. (Ha nil) can be derived from Eqs. (fT3TT6)) by 
exchanging Xi and X2- With ni = —1 and 71,2 = —1 
we have to flip the direction of xi and X2 ■ Finally, with 
ni = —1 and n2 = 1 'we have to exchange and flip the 
vectors. 

Suppose that the indices rii are randomly generated 
with probability pi{ni,n2) = 1/4. The joint probability 
distribution of Xi and Ui can be suitably written in the 
form 

p{xi,X2,ni,n2\b) = p{ni,n2\xi,X2;b)p^{xi,X2), (19) 



where 



p{ni,n2\xi,X2;b) = 9 (nib-yij 9 hi2b ■ y2J (20) 



and 



Pv{xi,X2) 



(47r)2 



(21) 



The model with probability distribution ()19|) and con- 
ditional probability ()18|) gives again the correct quan- 
tum probabilities, but now we have the nice property 
that the marginal probability distribution Pv{xi,X2) of 
the vectors Xi and X2 does not depend on the prepared 
quantum state b. This dependence is only in the con- 
ditional probability distribution p{ni, n2\xi,X2] b) of the 
discrete indices n^. This reformulation of the Kochen- 
Specker model gives a protocol with shared noise for 
simulating the communication of a qubit with just two 
bits of classical communication: Bob and Alice share the 
random vectors Xi and X2, generated according to the 
uniform probability distribution (l?T|) : given a quantum 

state &, Bob generates the discrete indices rii according 
to rule (I^Ul) : he sends them to Alice; given a measure- 
ment Ma, Alice generates the outcome a^ according to 
rule dUl). 



Any classical simulation of a quantum channel can be 
converted into a simulation of entanglement without in- 
creasing the amount of communication. Consider again 
the scenario of Sec. Ill Al with two qubits in the Bell 
state (H)). Alice and Bob perform two projective measure- 
ments. Ma = {01,02} and Mb = {61, &2}, respectively, 
and get outcomes Oq and bp. The marginal probabil- 
ity of /3 is uniformly distributed on the values 1 and 2. 
Furthermore, given Bob's outcome bp, Alice's outcome is 
generated as if she received the quantum state —bp di- 
rectly from Bob. Thus, the joint probability distribution 
of the outcomes can be simulated as follows. Bob ran- 
domly generates the outcome bp with uniform probability 
distribution p{f3) = 1/2. He then uses a classical model 
of a quantum channel for sending the quantum state —bp 
to Alice, who flnally generates her own outcome Cq,. The 
amount of communication of the derived entanglement 
model is equal to that of the quantum channel model. 

The classical model of a quantum channel previously 
derived from the KS model gives the following protocol 
for simulating entanglement. Alice and Bob share the 
random vectors, Xi and X2, and perform two projective 
measurements. Ma = {01,02} and Mb = {61,62}, re- 
spectively. Bob generates the outcome bp and the indices 
ni = ±1, ?i2 = ±1 according to the probability distribu- 
tion 

p{l3,ni,n2\xi,X2; Mb) = -9 (nibp ■ yA 9 f n2 6/3 ■ mj , 

He sends n^ to Alice. She generates the outcome 
Oq. according to the probability distribution deflned by 
Eq. (O. 

In this model the amount of communication is 2 bits. 
It is possible to reduce the communication cost by means 
of the transformation nij/i — >■ yi. Indeed the index ni 
becomes uncorrelated with the other variables and can be 
eliminated. Thus, we get a model with a communication 
cost of 1 bit, deflned by the conditional probabilities 

p{P,n\xi,X2]MB) ^9(bp- yi) 9 (nbp ■ ^2) , (22) 

p{a\xi,X2,n;MA) = 9[-aa ■ (yi + ny2)] . (23) 

That is. Bob generates the vector bp closest to j/i and 
sets the discrete index n equal to the sign of bp ■ y2- He 
sends n to Alice. She then generates the outcome Oq that 
is closest to yi +ny2. Note that this model is very similar 
to the Toner-Bacon model, with the only difference that 
the unit vectors Xi are replaced by the vectors yi, which 
are a linear combination of xt [see Eq. (jlSp ]. 

This model of entanglement can be put into a more 
synthetic form that will be useful in Sec. IIII Al Sup- 
pose that the measurement outcomes are Og and bs and 
the communicated index is n; then it is easy to show 



by Eq. 



that b 



'^ • 



> bp ■ Xn for any /3 e {1,2} 



and n £ {1, 2}. Furthermore from Eq. ((23)) we have that 
—Ua • Xn > —Sa ■ Xn for any a G {1,2}. Thus, the al- 
gorithm can be reformulated as follows. Bob evaluates 
the vectors bp € {61,62} and Xn G {a?i,a?2} that maxi- 
mize bp ■ Xn- He generates the outcome bp and sends the 
index n — {1,2} to Alice. Alice generates the outcome 
cia G {01,02} that maximizes —aa ■ Xn- This algorithm 
will be generalized to higher dimensions of the Hilbert 
space in the next section. 



III. GENERALIZING TO HIGHER 
DIMENSIONS 



In Sec. IIII Al the three approximate protocols for simu- 
lating entanglement in higher dimensions are introduced. 
In Sec. IIII Bj we will present a simple method for convert- 
ing an entanglement protocol into a protocol for simulat- 
ing quantum channels. 



A. Entanglement 

Let us consider the following scenario. Alice and Bob 
receive two A^-dimensional quantum systems in the en- 
tangled state 



IV'. 



AB, 



1 ^ 



mi 



(24) 



They perform local projective measurements using the 
set of orthogonal vectors { | t/fi ),..., |(/)7v)} = Ma and 
{IV'i), ..., I^at)} = Mb^ respectively. This scenario can 
be reproduced by local hidden variables augmented by 
some amount of communication. 

The gap between the classical and quantum algorithms 
in the HM problem implies that this communication can- 
not be smaller than e'^'-^^ in the case of bounded error, 
n — log2 N being the number of ebits. A stronger con- 
straint was given in Ref. |8|, where it was shown that 
0(2") bits of communication are necessary for an ex- 
act simulation. An exact two-way classical protocol with 
0(712") bits of communication on average was reported 
by Massar et al. [9]. Unlike the protocols considered 
in this paper, the protocol in Ref. [9!] does not require 
any local hidden variables. A one-way model of a quan- 
tum channel for a single qubit was recently reported in 
Ref. [155. It requires an infinite amount of communication 
for an exact simulation, but it has the nice property of 
encoding the communication in a single real variable, in- 
stead of two real variables, which are required for defining 
a quantum state. It can be easily converted into a clas- 
sical protocol of entanglement. Just as that in Ref. Q, 
this protocol does not need local hidden variables. 

In this section we present three approximate protocols 
for simulating entanglement that give accurate results for 
low values of N and require a communication of just n 



bits. They were numerically tested for a family of mea- 
surements. The maximal discrepancy between the best 
protocol and the quantum predictions is less than 1% for 
A^ = 3 and grows sublinearly with n{— log2 N) in the 
range numerically studied. The first protocol is a gener- 
alization of the Toner-Bacon model of a single Bell state, 
suitably readjusted in Sec. Ill Al The second protocol 
is a generalization of the alternative model presented in 
Sec. IIIBI The last model is similar to the second one, 
but with a different shared randomness. 
Protocol 1 

1. Bob and Alice share one random vector |a:) and a 
random basis {|yi), ..., |j/jv)} = Y . 

2. Bob generates the outcome |-0b) G Mb that is clos- 
est to \x). 

3. He sends Alice the index of the vector \ym) G Y 
that is closest to jf/'h). 

4. Alice generates the event |0a) G Ma that max- 
imizes {(l)*a\{\x){x\ + \yn){yn\}Wa). where |0:) = 

Note that this model gives the right marginal distri- 
butions for the outcomes of each party, since the vectors 
I a;) and \yk) are generated uniformly in the Hilbert space. 
This holds also for the other protocols. Thus, only the 
accuracy of the correlations has to be checked. 

Protocol 2a 



1 . Bob and Alice share a set R ; 
random vectors. 



{|a;i),...,|xAr)}ofA^ 



2. Bob evaluates the vectors \ipb) G Mb and |a::„) e R 
that maximize Ki/jblxn)!^. 

3. He generates the outcome 1-06) and sends n to Alice. 

4. Alice generates the event \(j)a) G Ma that maxi- 
mizes |(0*|x„)p. 

This protocol is the generalization of the classical 
model of entanglement derived in Sec. IIIBI As shown 
below, it is more accurate than protocol 1 for TV > 10, 
in the sense that the maximal discrepancy (in the set of 
tested measurements) is lower. This is true at least for 
A'^ < 32, which is the maximal dimension we considered 
in the numerical simulations. The last protocol is similar 
to protocol 2a, but the vectors \xk) are orthogonal. 

Protocol 26: the same as protocol 2a, but Bob and 
Alice share a random orthogonal basis. 

While this protocol is not exact for A^ = 2, it is more 
accurate than protocol 2a for A^ > 2. As shown by the 
numerical simulations, the accuracies of protocols 2a and 
26 approach each other in the limit of high dimensions. 

It is interesting to note that the three protocols share 
a common feature. Kochen and Specker proved that any 
deterministic hidden variable theory equivalent to quan- 
tum theory is contextual [lj|. Although our protocols 
are not completely equivalent to quantum theory, they 
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FIG. 1: Probability of the event |0i) given \'4)i) as a function 
of tp for A^ = 3. The solid line is the quantum prediction. 
The circles and crosses are the probabilities obtained from 
protocols 1 and 26, respectively. 



satisfy this general constraint, that is, the probability of 
a joint event (|(/)q,), |-0/3)) depends on the whole set of 
orthogonal vectors Ma and Mb- Indeed, all the vec- 
tors in the sets are involved in the maximization pro- 
cedures used by the protocols. It is interesting to note 
that a similar procedure of maximization was used in 
an approximate hidden- variable model of qutrit reported 
in Ref. pjj. The Kochen-Specker theorem was recently 
generalized to probabilistic theories in Ref. |17| . 

Since these protocols are contextual and approximate, 
we can expect that the probability distribution of two 
events \^p) and |0a) again depends on all the vectors in 
the sets Ma and Mb ■ For the sake of simplicity, here we 
will consider the set of one-parameter measurements 



l^i) = cos (^IV'i) - sin v3|i/'2) 
\(f)2) = sin (p|^i) + cos ip\'4'2) 
|0fc) = |Vfc)forfc = 3,...,7V, 



(25) 



that is, the vectors \(j)k) are set equal to IV'fc) for k > 2. 
The results do not change qualitatively with a different 
choice of the measurements. 

The quantum probabilities of the joint events \(j)b) and 
IV'a) for a 7^ 1,2 or b ^ 1,2 are obviously equal to Sa.b- 
The probabilities P{a, b\if) of the outcomes a = 1, 2 and 
6 = 1,2, given ip, are 



1 



P(l,l|^) = P(2,2|^) = -cosV (26) 

P(l,2|^)=P(2,l|^) = lsin2^. (27) 

Note that the measurement with 1^9 = i^sq is equiva- 
lent to the measurement with if = ^ — ifQ and |V'i) and 
IV'2) swapped. Furthermore, J2k i^i^'^l'fi) ~ ^- Thus, 
it is sufficient to evaluate the discrepancy between the 
model and quantum theory for outcomes a = b = 1 and 
ip — [0, 7r/2]. This discrepancy contains the full infor- 
mation about the discrepancy of any other event with 
constraint (ESl). 



In Fig. [1] we report the probability of an event |(/)i), 
given \ipi), as a function of 1^9 for A^ = 3. The solid line is 
the quantum prediction cos^ ip. The circles and crosses 
are the Monte Carlo results, generated by the protocols 
1 and 2b. The maximal discrepancy between protocol 
2b and cos^ (p is very small, less than 0.01, and is barely 
perceptible in the figure. Protocol 2a gives the worst 
results with a maximal discrepancy about 0.025. Its data 
are not reported in the figure. The maximal discrepancy 
for protocol 1 is about 0.014. The numerical simulations 
were performed with 10^ Monte Carlo executions, which 
make the statistical error inappreciable in the figure. 
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FIG. 2: Same as Fig. [l] with A^ = 4 (dashed line) and iV = 10 
(dotted line). 

The same simulations were executed for higher dimen- 
sions of the Hilbert space; the results for A^ = 4 (two 
qubits) and iV = 10 are plotted in Fig. [2j The proba- 
bilities generated by protocol 1(2&) are reported at the 
left-hand (right-hand) side. It is interesting to note that 
protocol 26 gives the right probability for (p — 0, tt/2, that 
is, when \(j>i) and |^i) are parallel or orthogonal. How- 
ever, a small discrepancy was noted for measurements 
that do not satisfy constraint (1^ . 

The maximum discrepancy of the three protocols is re- 
ported in Fig.[3]as a function of A''. The horizontal axis is 
in logarithmic scale. For large values of N, protocols 2a 
and 26 have the same discrepancy. This is because in pro- 
tocol 1 the vectors \xk), which are generated randomly, 
have a very high probability of being almost orthogonal 
to each other for TV ^ 1. This fact makes protocol 2a 
almost indistinguishable from protocol 26 in high dimen- 



B. Quantum channel 



In Ref. llj it was shown that any protocol that simu- 
lates entanglement of two qubits can be converted into a 
classical model of a quantum channel for a single qubit. 
A very simple method for converting a general classical 
model of entanglement into a classical model of a quan- 
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FIG. 3: Maximum discrepancy between the quantum predic- 
tion and the classical protocols as a function of A'^. The data 
are interpolated by lines, plotted as a guide of eyes. 



turn channel is as follows. In the entanglement model, 
Bob and Alice share a set R of random vectors. This 
means that they have a common list of noise realizations 
Rk, with k = 1,2,3,.... They start from k = 1 and 
at each execution of the Monte Carlo simulation they 
read the next element of the list. It is possible to con- 
vert the protocol into a protocol for simulating quantum 
channels by increasing the communication by n bits on 
average. Suppose that Bob receives the quantum state 
1^). He selects a measurement Mb={\iPi),---,\iPn)} so 
that j-^i) = \ijj)- Using the protocol for simulating entan- 
glement and the first noise realization in the shared list. 
Bob generates a vector \ipb}- If IV'ft) ¥^ \i^)^ Bob interrupts 
the protocol for simulating entanglement and reads the 
next realization of the noise in the shared list; he repeats 
the procedure until \ipb) = \ip)- He then executes the 
communication procedure as established by the entan- 
glement protocol. Furthermore he sends the number of 
noise realizations that Alice has to skip. This additional 
information is equal to log2 N = n on average. 

Thus, using this strategy, it is possible to convert the 
three models of entanglement into protocols for simulat- 
ing quantum channels. This conversion requires doubling 
the amount of communication on average. 



IV. CONCLUSION 

In this paper we have presented three approximate 
one-way communication protocols for simulating the out- 
comes of local measurements, performed on bipartite en- 
tangled states. These protocols use an amount of com- 
munication equal to the number n of ebits. We have 
seen that they can be converted into approximate pro- 
tocols for simulating quantum communication channels. 
Approximate models like these can be useful for detect- 
ing if a quantum communication algorithm can be effi- 
ciently simulated by some classical algorithm. This is 



the case when a quantum communication algorithm uses 
states and measurements that an approximate model of 
a quantum channel can efficiently simulate with zero or 
bounded error in any dimension. 

The results reported in this paper can be improved in 
different ways. It is interesting to note that the protocols 
2a and 26 have the same general structure, but different 
shared noises. The general structure is as follows. Bob 
and Alice share a set i? = {|xi), ..., |a;A')} of N random 
vectors with probability distribution p(\xi), ...,\xn))- 
Bob evaluates the vectors \ipb) £ Mb and |x„) £ R that 
maximize the function 



Fi{b,n) = \{iPb\xn)\' 



(28) 



He generates the outcome \i/jb) and sends n to Alice. Al- 
ice generates the event |(/)a) £ Ma that maximizes the 
function 



F2(a)^|(0:K)|^ 



(29) 



We have seen that a suitable choice of the noise distri- 
bution can considerably reduce the error. Indeed, in the 
three-dimensional case a change of noise dropped the er- 
ror from 2.5% (protocol 2a) to less than 1% (protocol 
26). 

Thus, protocols 2a and 26 can be improved by evaluat- 
ing the optimal probability distribution p(|xi),..., \xn)) 
that minimizes the nraximal or average error over a set of 
measurements. Of course, the error cannot be reduced to 
zero if the whole set of measurements is considered, since 
this would require an exponential amount of communi- 
cation. However, another possible improvement can be 
reached by augmenting the communication cost, namely 
log2 N. This is achieved by increasing the number N 
of random vectors \xk). We can expect that this strat- 
egy and the optimal choice of p will enhance the accu- 
racy of the model, as defined through the functions in 
Eqs. (|28l29p . For a sufficiently large amount of commu- 
nication, one can hope to make the error bounded. No- 
tice that, as said in the introduction, the HM problem 
establishes the lower bound 2(^^(^) for the communica- 
tion cost, whereas the best known protocols require an 
amount of communication that scales as n2" |6|, Q • ^t 
present it is not clear if such protocols are also optimal. 
This open problem can be investigated by studying the 
class of protocols introduced in this paper. In conclusion, 
our models can give some indication for finding optimal 
one-way communication protocols that classically simu- 
late quantum channels and entanglement. Furthermore, 
they can be used for testing the efficiency of a quantum 
communication protocol versus classical protocols. 
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